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ABSTRACT 

We develop a theory of buoyancy instabilities of the electron-ion plasma with the heat 
flux based on not the MHD equations, but using the multicomponent plasma approach. We 
investigate a geometry in which the background magnetic field, gravity, and stratification 
are directed along one axis. No simplifications usual for the MHD-approach in studying 
these instabilities are used. The background electron thermal flux and collisions between 
electrons and ions are included. We derive the simple dispersion relation, which shows that 
the thermal flux perturbation generally stabilizes an instability. There is a narrow region 
of the temperature gradient, where an instability is possible. This result contradicts to a 
conclusion obtained in the MHD-approach. We show that the reason of this contradiction is 
the simplified assumptions used in the MHD analysis of buoyancy instabilities and the role 
of the longitudinal electric field perturbation, which is not captured by the MHD equations. 
Our dispersion relation also shows that a medium with the electron thermal flux can be 
unstable, if the temperature gradient of ions and electrons have the opposite signs. The 
results obtained can be applied to ICM and clusters of galaxies. 
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1. INTRODUCTION 

Various instability mechanisms are studied to understand some of the main features and 
processes in the astrophysical objects depending on their physical properties. Convective or 
buoyancy instability arising as a result of stratification is among those instabilities that may 
operate under different circumstances from the stellar interiors (e.g., Schwarzschild 1958), 
accretion disks (Balbus 2000, 2001), and neutron stars (Chang & Quataert 2009) to the hot 
accretion flows (e.g., Narayan et al. 2000, 2002) and even galaxy clusters and intercluster 
medium (ICM) (e.g., Quataert 2008; Sharma et al. 2009; Ren et al. 2009). Analogous insta- 
bilities also exist in the neutral atmosphere of the Earth and ocean (Gossard & Hooke 1975; 
Pedlosky 1982). Diversity of the astrophysical objects, in which convective instabilities may 
have a significant role, leading to turbulence and anomalous energy and matter transport, is 
a good motivation to explore this instability either through linear analytical analysis or by 
direct numerical simulations from different physical point of views. Although the significant 
role of convection in the transport of energy in stellar interiors is a well-known physical 
process, theoretical efforts to understand convective energy transport in the tenuous and hot 
plasmas such as ICM (Sarazin 1988) have lead to some results over recent years. 

According to the standard Schwarzschild criterion, a thermally stratified fiuid is convec- 
tively unstable when the entropy increases in the direction of gravity (Schwarzschild 1958). 
By taking into account the anisotropic heat fiux in plasmas where the mean free path of ions 
and electrons is much larger than their Larmor radius, one obtains additional instabilities for 
short wave numbers with larger growth rates than that without thermal flux. These instabil- 
ities have been shown to arise when the temperature increases in the direction of gravity at 
the absence of the background thermal flux (the magnetothermal instability (MTI)) (Balbus 
2000, 2001) and when the temperature decreases along gravity at the presence of the latter 
(the heat buoyancy instability (HBI)) (Quataert 2008). Both MTI and HBI have been sim- 
ulated in 2D and 3D by many authors over recent years (e.g., Parrish et al. 2008; Parrish & 
Quataert 2008; Parrish et al. 2009). Following recent achievements in the convective theory, 
it has attracted attention of the authors for analyzing its possible role in ICM after a long 
time discounting. Majority of the mass of a cluster of galaxies is in the dark matter. How- 
ever, around 1/6 of its mass consists of a hot, magnetized, and low density plasma known 
as ICM. The electron density is n,. ~ 10^^ to 10^^ cm^'^ at the central parts of ICM. The 
electron temperature Tg is measured of the order of a few keV, though the ion temperature 
Tj has not yet been measured directly (e.g., Fabian et al. 2006; Sanders et al. 2010). The 
magnetic field strength B in ICM is estimated to be in the range 0.1-10 /xG depending on 
where the measurement is made (Carilli & Taylor 2002) which implies a dynamically weak 
magnetic field with p = STTUeTe/B^ ^ 200 - 2000. Thus, ICM with the ion Larmor radius 
10^~^ cm (Tj ~ Te) and the mean free path 10^^~^^ cm is classified as a weakly coUisional 
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plasma (Carilli & Taylor 2002). In simulating ICM, it is important to consider anisotropic 
viscosity as well because the Reynolds number is very low (Lyutikov 2007, 2008). Another 
important physical agent is cosmic rays. Recent studies show that centrally concentrated 
cosmic rays have a destabilizing effect on the convection in ICM (Chandran & Dennis 2006; 
Rasera and Chandran 2008). 

Theoretical models applied for study of buoyancy instabilities arc based on the ideal 
magnet ohydro dynamic (MHD) equations (Balbus 2000, 2001; Quataert 2008, Chang & 
Quatacrt 2009; Ren et al. 2009). Using of these equations permits us comparatively easily 
to consider different problems. However, the ideal MHD does not capture some important 
effects. One of the such effects is the nonzero longitudinal electric field perturbation along 
the background magnetic field. As we show here, the contribution of currents due to this 
small field to the dispersion relation can be of the same order of magnitude as that due 
to other electric field components. Besides, the MHD equations do not take into account 
the very existence of various charged and neutral species with different masses and elec- 
tric charges and their collisions between each others and therefore can not be applied to 
multicomponent systems. On the contrary, the plasma E-approach deals with dynamical 
equations for each species. From Faraday's and Ampere's laws one obtains equations for the 
electric field components. Such an approach allows us to follow the movement and changing 
of parameters of each species separately and obtain rigorous conditions of consideration and 
physical consequences in specific cases. This approach permits us to include various species 
of ions and dust grains having different charges and masses. In this way, streaming insta- 
bilities of rotating multicomponent objects (accretion disks, molecular clouds and so on) 
have been investigated by Nekrasov (e.g., 2008, 2009 a, 2009 b), which have growth rates 
much larger than that of the magnetorotational instability (Balbus 1991). In some cases, the 
standard methods used in MHD leads to conclusions that are different from those obtained 
by the method using the electric field perturbations. One of a such example is considered in 
Nekrasov (2009 c). 

In this paper, we apply a multicomponent approach to study buoyancy instabilities in 
magnetized electron-ion astrophysical plasmas with the background electron thermal flux. 
We include collisions between electrons and ions. However, we adopt here that cyclotron 
frequencies of species are much larger than their coUision frequencies. Such conditions are 
typical for ICM and galaxy clusters. In this case, as it is known, the heat flux is anisotropic 
and directed along the magnetic field hues (Braginskii 1965). We consider a geometry in 
which gravity, stratification, and the background magnetic field are all directed along one 
(z-) axis. In our approach, it is important to obtain exact expressions for species' velocities 
in an inhomogeneous medium. We give main equations and results. However, for those who 
are not interested in the mathematical details, they can directly refer to Sections 7 and 8. 
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The dispersion relation is obtained for cases, in which the background heat flux is absent or 
present. This gives a possibility to compare two cases. Solutions of the dispersion relation 
are discussed. 

The paper is organized as follows. In Section 2, the fundamental equations are given. 
An equilibrium state is considered in Section 3. Perturbed ion velocity, number density, and 
thermal pressure are obtained in Section 4. In Section 5, we consider the perturbed velocity 
and temperature for electrons. Components of the dielectric permeability tensor are found 
in Section 6. 



2. BASIC EQUATIONS 

We start with the following equations for ions: 

^= + gH E H Vj X B - Z/ie Vi - Ve , 1 

at miUi rrii rriiC 
the momentum equation. 



the continuity equation, and 



+ V • mvi = 0, (2) 



^ + • Vpi + 7p,V • Vi = 0, (3) 
ot 



the pressure equation. The corresponding equations for electrons are: 

= + QeE + — Ve X B - m^l^ei {^e " V^) , (4) 

rie c 

+ V • neVe = 0, (5) 
^ + Ve • Vpe + 7PeV ■ Ve = A - (7 - 1) V • Qe, (6) 

^+Ve-Vre + (7-l)7;V-Ve=--(7-l)-V-qe, (7) 

ot Tie Ue 

the temperature equation, where qe is the electron heat flux (Braginskii 1965). We neglect 
inertia of the electrons. In Equations (l)-(7), qj and rrij are the charge and mass of species 
j — i,e, Vj is the hydrodynamic velocity, rij is the number density, pj = njTj is the thermal 
pressure, 7} is the temperature, i/ie (lyei) is the collision frequency of ions (electrons) with 
electrons (ions), g is gravity, E and B are the electric and magnetic fields, c is the speed of 
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light in vacuum, and 7 is the adiabatic constant. We assume the electrons to be magnetized 
when their cyclotron frequency oUce = qeB/nieC ^ Ugg, where z/ge is the the electron-electron 
collision frequency. In this case, the electron thermal flux is mainly directed along the 
magnetic field, 

qe = -Xeb(b-V)Te, (8) 

where Xe is the electron thermal conductivity coefficient and b = B/i? is the unit vector 
along the magnetic field (Braginskii 1965). The term A compensates the temperature change 
as a result of the equilibrium heat fiux. Wc take only into account the electron thermal 
conductivity by equation (8), because the corresponding ion conductivity is considerably 
smaller (Braginskii 1965). 

Electromagnetic equations are Faraday's law 
and Ampere's law 

V X B (10) 
c 

where j — '^j QjnjVj. We consider the wave processes with typical time-scales much larger 
than the time the light spends to cover the wavelength of perturbations. In this case, one 
can neglect the displacement current in Equation (10) that results in quasineutrality both in 
electromagnetic and purely electrostatic perturbations. The magnetic field B includes the 
background magnetic field Bq, the magnetic field Bocur of the background current (when it 
presents), and the perturbed magnetic field. 



3. EQUILIBRIUM STATE 

At first, we consider an equihbrium state. We assume that background velocities are 
absent. In this paper, we study configuration, in which the background magnetic field, 
gravity, and stratification are directed along the 2;-axis. Let, for definiteness, g be g = —zg, 
where g > and z is the unit vector along the z-direction. Then, Equations (1) and (4) give 

1 dpio Qi , . 

Qi ^ -K- ^9 Eo, 11) 

miTiio oz rrii 

Qe — ^ — -C/O, (l^j 

miUeo oz rrii 

where (and below) the index denotes equilibrium values. Here and below we assume that 
qi — —qe- We see that equilibrium distributions of ions and electrons influence each other 
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through the background electric field Eq. In the case njo = rigo ( this equality is satisfied 
for the two-component plasma) and T^o — T^q, we obtain gi — = g/2. Thus, we have 
£^0 = nT'id/'^Qi- The presence of the third component, for example, of the cold dust grains 
with the charge and mass ^ rrii results in other value of Eq = rrifig/qd- In this 
case, the ions and electrons are in equilibrium under the action of the thermal pressure and 
equilibrium electric field, being gi ~ —gg. 



where 



4. LINEAR ION PERTURBATIONS 

Let us write Equations (l)-(3) for ions in the linear approximation, 

^— = 1 h Fa H Vii X Bo, (13) 

dnn driio , 

+ Viu-^ + riioV ■ Vii = 0, (14) 

-^ + 'yiu-^+7PioV-Va = 0, (15) 

Fa = —El - i^ie {vii - Vei) , (16) 



and the index 1 denotes the perturbed variables. Below, we solve these equations to find the 
perturbed velocity of ions in an inhomogeneous medium. 



4.1. Perturbed velocity of ions 

Applying the operator d/dt to Equation (13) and using Equations (14) and (15), we 
obtain 

^= -9iVva. + ^ [(7 - 1) (Vp.o) + W^oV] V ■ + ^ + x Bo. (17) 

at^ miTiiQ at rriiC dt 

We can find solutions for the components of Vji. For simplicity, we assume that d/dx = 0, 
because a system is symmetric in the transverse direction relative to the z-axis. The x- 
component of Equation (17) has the form 

—^=Filx + l^ciVilyi (18) 
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where oud — qiBo/rriiC is the ion cyclotron frequency. For the y-component of Equation (17), 
we obtain: 



-9i 



+ 4— V-Va + 



dt 



dt 



^2 dy ^^dy 

Here, Cgi = {'yTio/rrii)^^'^ is the ion sound velocity. 

Using Equation (18), a relation for Vny is given from Equation (19) as follows 

Then from Equation (18), we obtain 



(19) 



d 



dy 



(20) 



(21) 



Here, the following notations are introduced: 

Pii = -QiViiz + 4V • Va, 



Q 



ilx 



ci-'^ ily 



dt ' 



(22) 
(23) 
(24) 



The value Pn defines the pressure perturbation (Eq. [15]). We see from Equation (21) 
that when d/dt <^ the thermal pressure effect on the velocity Vnx is much larger than 
that on Viiy. The 2;-component of Equation (17) takes the form 



d ( dviu 



dtXdt 



Fiu 



dviiz , 
dz 



„2 d 



(1 -7)5i + c,^— 



(25) 



Let us now find V • Vji through vn^- Differentiating Equation (20) with respect to y 
and using expression (22), we obtain 



where the following operators are introduced: 



dQ 



ily 



dy 



(9^ d'^ 

J _ ^ I 2 _ 2 ^ 

' ~ dt^ " ""''dy^' 



d^ 



d 



-9i 



(26) 

(27) 
(28) 
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We can derive an equation for the longitudinal velocity vnz, substituting V • Vji found 
from Equation (26) in Equation (25), 

LsVa.=L^^ + L,^, (29) 

where operators L3 and L4 have the form 

- ^ J ~ w ^d^yd^ '''''''''' ^ ^ 

fd^ 2^ 5 2 dL. ,2 
+^9i [^+^ci]^ + Csi-r^M + (1 - 7) 9i 



,..(1_,),. + 4(|_||). ,31) 

For obtaining expression (30), we have used expressions (27) and (28). 

It is easy to see that at the absence of the background magnetic field and without taking 
into account electromagnetic perturbations (the right hand-side of Eq. [29]), the equation 
LsViiz — describes the ion sound and internal gravity waves. In this case, a sum of the last 
two terms on the right hand-side of expression (30) is equal to — c^ja;^j^, where is the 
(ion) Brunt- Vaisala frequency equal to 

(32) 

However, we see the existence of the background magnetic field considerably modifies the 
operator L3. Note the right hand-side of Equation (29) describes a connection between ions 
and electrons through the electric field Ei and coUisions. 



9i_ 
^2 



4.2. Specific case for ions 

So far, we have not made any simplifications and all the equations and the expressions 
are given in their general forms. Now, we consider further perturbations with a frequency 
much lower than the ion cyclotron frequency and the transverse wavelengths much larger 
than the ion Larmor radius. Such conditions are typical for the astrophysical plasmas. 
Besides, we investigate a part of the frequency spectrum in the region lower than the ion 
sound frequency. Thus, we set 

2 92 2 2 
^ci > c, — ; c,— > — . (33) 
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In this case, operators (27), (28), (30), and (31) take the form 



(34) 



d 



d 



(l-7)ft + c2 



d_ 

dz' 



Also, the operator L3 can be written for a case in which 

92 



(35) 



The small corrections in operators L3 and L4 are needed to be kept because some main terms 
in expressions for ion and electron velocities are equal each other (see below). Therefore, 
when calculating the electric current these main terms will be canceled and small corrections 
to velocities will only contribute to the current. 

For the cases represented by inequalities (33) and (35) when the operators have the form 
(34), the equations for and V • Vji become 



d 



""''dz dz dt^ 



d 92 



Bt 

Bz ul^By 



Bz 



dQ 



ily 



(36) 
(37) 



4.3. Ion velocity in the Fourier transform 

Calculations show that some main terms in expressions for vn^ (when calculating the 
current), V • v^i and Pn are canceled. Therefore, the small terms proportional to inhomo- 
geneity must be taken into account. To make this correctly, we can not make the Fourier 
transformation in Equations (36) and (37) to find the perturbed ion pressure Pn. However, 
firstly, we should apply the operator B/Bz to this variable for using Equation (36). It is 
analogous to obtain the term Bc^/Bz in expression (32) for the Brunt- Vaisala frequency. 
After that, we can apply in a local approximation the Fourier transformation assuming the 
linear perturbations to be proportional to exp(ikr— icut). As a result, we obtain for the 
Fourier-components vnzk, k • Vnk and Puk, where k — (k,^), the following expressions: 
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UJ 



1-i 



19i 



Fiizk + i 



Q 



ilyki 



(39) 



Pi 



Ik 



^ TP ■ ^ 

-rPilzk - '^yTa 



ilzk 



(40) 



+i 



^z^%^ci 



del 



(7-1)5^ + ^^-^- 



dz 



9i 



Qilyk- 



In expressions (38) and (39), we have omitted additional small terms at Qnyk, which 
are needed for calculation of Puk- When calculating the current along the 2;- axis, the main 
term ~ Qnyk in Equation (38) will be canceled. The contribution of the first term ~ Fn^k 
to this current has, as we shall show below, the same order of magnitude for the buoyancy 
instabilities as that of the term ~ QiQuyk- The same relates to expressions (39) and (40). 
Thus, the longitudinal electric field perturbations must be taken into account. However, in 
the ideal MHD, this field is absent. We see from expressions (38) and (39) that V ■ Vji ~ 
{gi/clj)viiz. This relation is the same as that for the internal gravity waves in the Earth's 
atmosphere (e.g., Nekrasov 1994). Using expression (40), we obtain velocities vnyk and vnxk 
from Equations (20) and (21), correspondingly. 



4.4. Perturbed ion number density and pressure 



It is followed from above that V ■ Va ^ 0. Let us find the perturbed ion number density 
and pressure in the Fourier- representation. From Equations (14), (38) and (39), we obtain 



nuk 



i J 2 Piizk i 
l^zCgi 



(7-1)9. + ^ 



Q 



ilyk- 



(41) 



Equation (15) gives dpa/dt — —minioPn. Thus, we obtain, using Equation (40), 

del 



Pilk • 7 p , lky9'. 

— —^ilzk + 



PiO 



kzCii^^ci 



dz 



UJ 



9i 



lyk- 



(42) 



Comparing Equations (41) and (42), we see that the relative perturbation of the pressure 
due to the transverse electric force Qnyk is much smaller than the relative perturbation of the 
number density. However, these perturbations as a result of the action of the longitudinal 
electric force Fn^k have the same order of magnitude. Thus, pnk/Pio ~ nuk/riio. This result 
contradicts a supposition pnk/Pio ^ nak/riio adopted in the MHD analysis of buoyancy 
instabilities (Balbus 2000, 2001; Quataert 2008) because the latter does not take into account 
the longitudinal electric field perturbations. From the results given below, it is followed that. 



as we have already noted above, the both terms on the right hand-side of Equation (41) have 
the same order of magnitude. 



5. LINEAR ELECTRON PERTURBATIONS 

Equations for the electrons in the linear approximation are the following: 

= \ h Fei H Vei X Bo, (43) 

UeO UeO UgQ C 

+ Velz + rieO V ■ Vel = 0, (44) 

+ ^^^^-Q^ + 7PeoV ■ Vei = - (7 - 1) V • q^i, (45) 

^ + ^eU^ + (7 - 1) TeoV • Vei = - (7 - 1) —V • q,i, (46) 

dTeO , dT^l , 9TeO 

Qel = - "oXeO-Q^ - OoXel-Q^, (47) 

Fei = 5eEi - rUel^ei (^el " Vji) ■ (48) 

Here, Xei — 5XeoTei/2Teo (and Xe ~ Te^^, see Spitzer (1962)) is the perturbation of the 
thermal flux conductivity coefficient. The perturbation of the unit magnetic vector bi is 
equal to bi^^y — Bi^ y/BQ and hiz — 0. The thermal flux in equilibrium is Qeo = — boXeo^^- 

We have seen above at consideration of the ion perturbations that the terms ~ l/i?^, 
where H is the typical scale height, are needed to be kept (see the last term in Equation 
(40)). Therefore, these terms are kept also for the electrons. 



5.1. Equation for the electron temperature perturbation 

Let us find equation for the electron temperature perturbation. The expression V ■ qg;^, 
where qei is defined by (47), is given by 

^ ^ dqely ^ dqelz _ dT^Q 1 dBiy d'^Tei ^ dXepdT^l d'^XeO rjn /^gx 

'^^ dy dz "^^^ dz Bq dy dz"^ dz dz dz'^ 
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Substituting this expression into Equation (46), we obtain 

DlTei = -Velz^ (7-1) ^eoV ■ Vgl + (7 - 1) ^^"^^ (50) 

oz Ueo oz Body 



where the operator Di is defined as 



- (7 - 1) XeO^ + 2^—^ + 



dt Ueo V dz^ dz dz dz'^ 



(51) 



5.2. Perturbed velocity and temperature of electrons 

We find now equations for components of the perturbed velocity of electrons. The 
x-component of Equation (43) has a simple form, i.e. 

Vely = F^lx, (52) 

rrieUJce 

where cuce = qeBo/meC. Applying the operator d/dt to the ^/-component of Equation (43) 
and using Equations (45) and (49), we obtain 



a, . Velx —Fely ] = " (T " 1) — p (53) 

where 



dt J dy ^ 



Pel = -QeVelz + C^^V ■ Vgl (54) 

and c^g = 7Peo/ fnineo- The variable Pel is analogous to Pa (see Eq. [22]), which defines the 
ion pressure perturbation. But for electrons, their pressure perturbation is also affected by 
the thermal conductivity (see Eq. [45]). 

Let us express V • Vgi through Veu, using Equation (52), 

dVelz 1 9Felx frr\ 

V ■ Vel = 55) 

dz rUeUce dy 
The ;2-component of Equation (43) takes the form 

0= ^ \ ^ h Feu. (56) 
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We consider further perturbations with the dynamic frequency d/dt satisfying the fol- 
lowing conditions: 



Ueo Oz'^ ^ dt ^ Ueo dz dz 



(57) 



In this case, the terms proportional to dxeo/dz in the temperature equation (50) (see 
[51]) are unimportant because the necessary small corrections proportional to d/dt in this 
equation will be larger than that ~ dxeo/dz. Thus, an inhomogeneity of the thermal flux 
conductivity coeflicient and its perturbation can be neglected. We further apply the operator 
d/dt to Equation (56) and use Equations (44), (45), (49), and (55). As a result, we obtain 



d_ 

'Wz 



dv, 



elz 



dz 



_dF^ 
rriidt 

+ (7-1 



+ 



XeO 



1 - 7) 5e + cl 



d_ 

dz 
1 d^B.y 



1 dFeu 
dy 

d'T,, 



(58) 



+ 



^j^eo \ dz Bq dydz dz^ 



Equation for the temperature perturbation under conditions (57) has the form 



(7-1: 



Xeod^ 

rieo dz^^ 



d_ 
dt 



el 



Velz^TT^ + (7 - l)^eO 



dz 
-(7-1) 



dVeU 

dz 



1 dF, 



elx 



dy 



(59) 



XeO dTeo dBiy 



Ueo dz Body' 

where we have used Equation (55). Substituting Tgi in Equation (58) and carrying out some 
transformations, we find equation for the longitudinal velocity Veu 

d^eiz S^Fei, Ueo fd_y'd^F,u . 1 d^F,,, 

dz) 

1 



dz^ 



Teodzdt XeO 



+— 

'^se 



del 



19e 



Teodt'^ meOJce dydz^ 

elx 



(60) 



d^F.: 



dz I nieUJce dydz T^odz Bq dydt 

The correction proportional to dF^ix/dt is absent. The last term on the right hand-side of 
Equation (60) is connected with the background electron thermal flux (Quataert 2008). 

Prom Equations (59) and (60), we can find equation for the temperature perturbation 



/ -.^XeO d 

7-1 — — 
rieo dz 



d^T,, ^ dTeodB^y 
dz^ dz Body 



eO 



(7 

-(7-1) 
dT,o 



d(? 



dFei 



1 dF, 



elx 



rrieUJce dy 



(61) 



£ _ "[ho 

dt ^ XeO 



d_ 
dz 



dt^ 



-7- 



dz 



d_ 
dz 



-1 



d^B 



Bodydf 
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It is followed from results obtained below that all terms on the right-hand side of Equation 

(61) (except the correction ~ d'^Feiz/dt^) have the same order of magnitude (see Section 
[4.3]). The left-hand side of this equation is larger (see conditions [57]). Thus, the temper- 
ature perturbation in the zero order of magnitude can be found by equaling the left part of 
Equation (61) to zero. However, the right side is necessary for finding the transverse velocity 
perturbation Veix 

To find the velocity feiz, we need to calculate the value P^i (see Eqs. [53] and [54]). 
Performing calculations in the same way as that for ions (see Section 4.3), we obtain 



Oc^ 

cL^ + (7-l).e + ^ 



dFeU ^ dVel 

dz 



+9e 



(7-1)...^ 



rriidt 

dFelx 



dy 



(62) 



where we have introduced the notation connected with the thermal flux. 



Ki = (7-1) 



XeO pTeO 1 dB^y ^ g^Tei 

miUeo \ dz Bq dy dz"^ 



(63) 



Equation (62) can be re-written in the form, which is convenient for finding the velocity 
Veix- Using Equation (61), we obtain 



dz' 



{Pel - Kl) 



7 



rriidzdt c^g 



dc^ 

{l-l)9e+ - 



dz 



dF^ 
m.}dt 



dc^ 

(7-1).. + ^ 



19e + 



dd 



1 dF, 



elx 



dz J meUJce dy 



dd 



dT,. ( d 



-1 i^i 



d^B 



Teodz \dz J Bodydt 



(64) 



It is easy to see that Equation (53) has the form 

1 „ 



d 

Ft ' "^^^ " 



rrieiOce 



ely 



1 d 

^(Pel-Kl) 



(65) 



Thus, the main contribution of the fiux described by Equation (63) does not infiuence on 
the electron dynamics. Applying to Equation (65) the operator d"^ jdz^ and using Equation 
(64), we find an equation for the velocity v^ix- 
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6. FOURIER CURRENT COMPONENTS 



6.1. Fourier velocity components of ions and electrons 



Let us give velocities of ions and electrons in the Fourier-representation. From Equations 
(20), (21), and (40), we have 



1 



CO' 



1 + 



1 H T Qiixk + i-p; 



1 9 3 ^ilyk~ 



1 ky (^_.a^^^^ 



U!r.i kz 



/t2 



Izk-i 



k^. 



, U ky f ■'^■i \ TP 

Qilyk + ^~2"ir [ ^ ~ ^JT I ^i^zk- 



Here and below, we have introduced notations 

1 



"he — 2 

si.e 



(7 - 1) gi,e + 



dc 



2 

si,e 



dz 



The velocity Vnzk is given by Equation (38). 
From Equations (64) and (65), we find 



'^elxk — ^ 



.tteCtkl 



uujci kl V m,A^,,, 



Op. ^elxk + UJ , 5— 

k^T^odz Bq 

1 



+ Feiyk-^ 1 - ^7T^ / 

rUeUJce kz \ kzj nieUJce 



■Felzk -I 



where the following notation is introduced: 

1 



del 



2 



Equation (60) also gives us 



Velzk — 



ky 1 _ 



, - elxk " , 9 1 

kzUleUJce k^ \ nieUJce 



.ky ( 1 9Teo B\yk 

KT^^dz Bo 



—I 



klTeo 



1 + l''^ Felzk- 

XeQl^z J 



The velocity v^xy is defined by Equation (52). 



(66) 
(67) 



(68) 



(69) 



(70) 



(71) 
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6.2. Fourier electron velocity components at the absence of heat flux 

To elucidate the role of the electron thermal flux, we also consider the dispersion re- 
lation when the flux is absent. Therefore, we give here the corresponding electron velocity 
components: 

Velxk — ~^T9 ^elxk-\ ^ elyk ~ T \ ^ ~ ^ elzki \' ^) 

k^UJUJcimeUJce rUeljJce \ k^ J TTleUJce 

Velzk^T~ 1 ~ T Pelxk-ijlT^ 1 " i-, J- ^elzk- (73) 

Comparing expressions (69) and (71) with these equations, we see that the thermal flux 
under conditions (57) essentially modifies the small terms in the electron velocity. 



6.3. Fourier components of the current 

We find now the Fourier components of the linear current ji = qiUioVn + qeUeoVei- It 
is convenient to consider the value At^i^i/u. Using expressions (38), (52), and (66)-(71), we 
obtain the following current components: 

jlxk — (IxxEixk + idxyEiyk — dxzEizk (74) 



~bxx (^ilxfc ^elxfe) ibxy {Vilyk ^elyfe) ~l~ ^xz {t^ilzk ^elzfe) ) 



47rz 

'jlyk — ~iO'yxElxk + dyyEiyk — ttyzEizk (75) 



-\-ibyx {Vilxk ^elxfc) byy (fjli^fc l^elyk) ~l~ byz (l^ilzfc l^elzk^ ) 



47ri 

-jlzk — —ClzxEixk ~ ClzyEiyk + dzzEiz (76) 



+bzx (Vilx - Velx) + bzy (vny - Vely) - b^z (Vilz " Velz) 



When obtaining expressions (74)-(76), we have used notations (16), (23), (24), and 
(48) and equalities qe = —qi, Ueo — riio, meVei — triiPie- We have also substituted Biy^ by 
{kzc/uj)Eixk (see below). The following notations are introduced above: 

fe' f. kl Qjaj + a,&,cg, kl g.cg, dT*^ \ 
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and 



_ _ _^ ^ 9i(^i \ _ ^pi ky ^ 



^yx 






Cl 


U2 




1 ' 






^pi ky 

^li kz 






1 




U2 






) 



UUcA kl 



bp'- h f n- Fit* 

— _Zi n — n — _t^'JJL n — y \ h — _ -I- 

ayy - , ,2 . ^ry ^ ,2 ,^ , ^2 I '^^ ^ T^n^^ 



SI 



^li^ie rrij k^ / _ k^gia4+_aecj,^be 
^ci 1i ^z V k'^ 



= -^-^ r-"M , (78) 



, _ ^pi fey / _ 5fi_\ m£ 
h - 

Here Upi = {Annioqf/mi) is the plasma frequency and k"^ — ky+k^. The terms proportional 
to T*Q are connected with the background electron thermal flux. 

Calculations show that to obtain expressions for aij without thermal flux, using electron 
velocities (72) and (73), we must change be by ge/clg, put T*q = 0, and take 7 = 1 in terms 
a^z and a^^. 



6.4. Simplification of collision contribution 

From the formal point of view, an assumption that electrons are magnetized has only 
been involved in neglecting the transverse electron thermal flux. In other respects, a rela- 
tionship between cuce and i/gi or cUd and 1/^^ (that is the same) can be arbitrary in Equations 
(74)-(76). We further proceed by assuming that u <^ Ud. In this case, we can neglect the 
coUisional terms proportional to b^y and by^- However, the system of Equations (74)-(76) 
stays sufficiently complex to find ji through Ei. Therefore, we further consider the specific 
case in which the frequency uj and wave numbers satisfy the following conditions: 



4 k'^ ^ie 

where 

2 2 

2 ^si^se 



vik^^ v,e^ klH^Ll' ^^^^ 
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It is clear that conditions (79) can easily be realized. In this case, the current components 
are equal to 

47ri 

jlxk ^xx-^lxk ~l~ ^^xy-^lyk ^xz-^lzki (^-'■) 

Airi . • 771 77 771 

Jlyk — '^^yx-t^lxk ~l~ ^yy^lyk ^yz^lzki 

UJ 

Am , 

Jlzfe — ^zx-^lxk ^zy-^lyk ~l~ ^zz-t^lz- 

UJ 

Components of the dielectric permeability tensor eij are the following: 

_ ■ i^ie ky {aj - -fae) _ Vie ky {aj - 'jae) , . 

Sxx - a^x + ^2 _ .^^^ a,,, e,y - a^y + (i _ ^^^) 

o-j.^ ujUi^ ky (^zx 



cjikAi-idzY 

Ojyz ^zx Ojzy 



~ {l-^dzy''' 
where we have used notations (78) 

Parameter dz defines the coUisionless, ^ 1, and coUisional, ^ 1 regimes. Below, we 
derive the dispersion relation. 



7. DISPERSION RELATION 

Prom Equations (9) and (10) in the Pourier-representation and using system of equations 
(81), we obtain the following equations for the electric field components: 

{v? — Exx) Eixk — iSxyEiyk + SxzEizk = 0, (84) 
i^yxEixk + {pi — Syy) Eiyk + {—nyUz + Syz) Eizk = 0, 
^zxE\xk + {—l^yl^z + ^zy) Eiyk + {riy — Szz) Ei^k — 0, 

where n = kc/u;. The dispersion relation can be found by setting the determinant of the 
system (84) equal to zero. In our case, the terms proportional to Sxy and Syx can be neglected. 
As a result, we have 

{jl ^xx) \j^y^yy ~l~ ij^z ^2/1/) ^-z-z f^yf^z i^yz ~l~ ^zy^ ~l~ ^yz^zy^ 



(85) 



The above dispersion relation can be studied for different cases. In subsequent sections, 
we consider botli tlie collisionless and collisional cases. 



7.1. Collisionless case 



We assume now that the condition 



«1, 



(86) 



is satisfied. Then, using notations (77) and (82), the dispersion relation (85) becomes 



where ca — Bq/ {A-nrriiniQY/'^ is the Alfven velocity and 



(87) 



For obtaining Equation (87), we have used the condition kyc'l/uj'^- <C f. We see that 
there are two wave modes. The first wave mode, = k'^c\, is the Alfven wave with a 
polarization of the electric field mainly along the y-axis (the wave vector k is situated in the 
y — z plane) . This wave does not feel the inhomogeneity of the medium. The second wave 
has a polarization of the magnetosonic wave, i.e. its electric field is directed mainly along 
the X-axis (see below). This wave is undergone by the action of the medium inhomogeneity 
effect. The corresponding dispersion relation is 



u^^klcl + n^'-;^. 



(89) 



The expression (88) can be further simplified using equations (11), (12), (68), (70), and 
(80). As a result, we obtain 



Q2 



7 



(7 - 1) rriig + 7 



dz 



mg + 



dz 



(90) 



We have pointed out at the end of Section (6.3) what changes must be done in expressions 
(77) and (78) to consider the case without heat flux. This case follows from Equation (90), 
if we omit the term d (Tgo + T*q) / dz and put 7 = 1 in the first multiplier. Then becomes 
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^si ~ "-'se 



(7-1)^ + 



seJ 



dz 



(91) 



This is the Brunt- Vaisala frequency. Comparing (90) and (91), we see that the heat flux 
stabihzes the unstable stratification. The presence of the background heat flux does not play 
of principle role. If the temperature decreases in the direction of gravity {dTi^eo/dz > 0), a 
medium is stable. Solution (90) describes an instability regime only when 

7-1 9^0 1 

27 ^ dz 2 ^ 

where (Tjo ~ Tgo — To). We also note that Jl^ can be negative if gradients of Tjo and Teo 
have different signs. 

For a comparison, we give here the corresponding dispersion relation by Quataert (2008) 

'd\nTo\ kl 



2 



dz J k'^^ 



which is discussed in Section 8. 



7.2. Collisional case 



We proceed with the collisional case when 



> 1. (92) 



k^cl 

In this limiting case, we obtain again Equation (89) 



7.3. Polarization of perturbations 

Let us neglect in the system of equations (84) the small contributions given by Exy and 
Syx- Then, for example, in the coUisionless case, we obtain for the second wave u'^ ^ k^c^, 

k 

E\yk — -j-Eizki (93) 

Eizk — —Eixk ^ Eixk- 



Thus, the second wave has a polarization of the electric fleld mainly along the x-axis. 
In spite of that the component Eizk ^ Ei^ki it is multiplied by a large coefficient in the first 
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equation of the system (84). As a result, the contribution of this term is the same on the 
order of magnitude as that of the first term. 

In the coUisional case, the component Ei^k is also defined by Equation (93). However, 
its contribution to the first equation of the system (84) can be neglected. 



8. DISCUSSION 

Dispersion relation (87) with defined by Equations (88) or (90) considerably differs 
from that given in (Quataert 2008) for the case of our geometry. The reason goes back to 
the assumptions made in the MHD analysis of buoyancy instabilities pi/po <^ Pi/ Po, where 
p and p denote the pressure and mass density of fluid, and the condition of incompressibility 
V • Vi = 0, where Vi = Vji is the perturbed fluid velocity. We now shortly show how one 
can obtain the result of Quataert (2008) in our geometry, using these assumptions. We sum 
Equations (13) and (43) and use the Ampere's law (10). The components of the equations 
become. 



dt 4:71 Po dz 

dviiy _ dpi Bq f dBu dBiy 



(94) 



dt pody 47rpo V dy dz 

dvju ^ dpi _ nei 
dt Pody no ' 

where riio = Ueo = no, po = niino, and we have used nn = nei- The components of the ideal 
magnetic induction equation are the following: 

dBix ^ dviix /„^s 

dBiy ^ ^ dvjiy 

dt ^ dz ^ 

dBi^ _ _^ dvjiy 

dt ^ dy 

We note that the last equation (95) is a consequence of the second equation (95) and 
V • Bi = 0. 

The first equations of the systems of equations (94) and (95) describe the Alfven waves, 
which are split from the other perturbations (see Eq. [87]). Applying operators d^/dz^dt 
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and d^/dydzdt to the second and third equations of the systems of equations (94), corre- 
spondingly, using equation V • Vji = and the second equation (95), and subtracting one 
equation from another, we obtain 



^ dy'^ dz"^ J dt"^ \ dy^ dz^ J dz^ uq dydzdt 
Also, from Equation (61), we have (see conditions [57]), 



dz'^ dz Body 



(97) 



Taking into account that T^i/TgQ = — ngi/ngo, differentiating Equation (97) over t, using 
the second equation (95), and substituting the equation obtained in Equation (96), we find 



By neglecting the contribution of the magnetic field, this equation coincides in the Fourier- 
representation with Equation (22) in Quataert (2008). 

However, the presence of the longitudinal electric field perturbations Ei^ results in that 
Pi,eik/Pi,eo ~ ''ii,eife/'«'i,eo (for ions, scc Scction [4.4]). Besides, Equation (97) together with 
equation V • Vji = lead to the nonphysical equation 



Vel2 



eO 



dt dz 
Thus, Equation (98) is incorrect. 

From the dispersion relation (85), we see the necessity of involving the contribution of 
values e^z, £zx, and e^z in the coUisionless case (86) (values e^z and e^^ give the last term 
on the right hand-side of Eq. [88]). This means that contribution of currents jix ~ -Biz and 
jiz ~ Eix, must be taken into account. However, the role of the longitudinal electric 
field £"12 in the MHD equations is not clear. The same also relates to the coUisional case 
(92). In the current jixk-, we must take into consideration the contribution of the current 
jizk as a result of collisions, which is proportional to Ei^k (see Eqs. [74] and [76]). 

Thus, the standard MHD equations with simplified assumptions are not applicable for 
the correct theory of buoyancy instabilities. Such a theory can only be given by the multi- 
component approach used in this paper. 

The results following from Equation (90) show that the thermal fiux stabilizes the 
buoyancy instability. The instability is only possible in the narrow region of the temperature 
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gradient (see Section [7.1]). The presence of the background electron thermal heat (the term 
~ T*q) does not play an essential role. An instability is also possibly, if the temperature 
gradients of ions and electrons have the opposite signs. 

The contribution of collisions between electrons and ions depends on the parameter dz 
defined by Equation (83). In the both limits (86) {d^ <^ 1) and (92) {dz » 1), the dispersion 
relation has the same form. 

We would like to say a few words about the Schwarzschild criterion of the buoyancy 
instability. It is generally accepted that this instability is possible, if the entropy increases 
in the direction of gravity. From a formal point of view, it is correct, if we take the Brunt- 
Vaisala frequency N in the form (e.g. Balbus 2000), 

^2 ^ 1 dp dlnpp-^ 

7P dz dz 

However, this expression can easily be transformed to expression (32). Thus, we see that for 
instability to exist, the temperature must increase along gravity and exceed the threshold. 



9. CONCLUSION 

In this paper, we have investigated buoyancy instabilities in magnetized electron-ion 
astrophysical plasmas with the background electron thermal flux, using the E-approach when 
dynamical equations for the ions and the electrons are solved separately via electric field 
perturbations. We have included the background electron heat fiux and collisions between 
electrons and ions. The important role of the longitudinal electric field perturbations, which 
are not captured by the MHD equations, has been shown. We showed that the previous 
MHD result for the growth rate in the geometry considered in this paper when all background 
quantities are directed along the one axis is incorrect. The reason of this has been shown to 
be in simplified assumptions made in the MHD analysis of the buoyancy instabilities. 

We have adopted that cyclotron frequencies of species are much larger than their coUision 
frequencies that is typical for ICM and galaxy clusters. The dispersion relation obtained 
shows that the anisotropic electron heat flux, including the background one, stabilizes the 
unstable stratiflcation except the narrow region of the temperature gradient. However, when 
gradients of the ion and electron temperatures have opposite signs, the medium becomes 
unstable. 

Results obtained in this paper are applicable to the magnetized weakly coUisional strat- 
ified objects and can be useful for searching sources of turbulent transport of energy and 
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matter. It has been suggested that buoyancy instabihty can act as a driving mechanism to 
generate turbulence in ICM and this extra source of the heating may help to resolve cooling 
flow problem. However, all previous analytical or numerical studies are restricted to the 
MHD approach. Our study shows that when the true multifluid nature of the system with 
the electron heat flux is considered, one can not expect buoyancy instability unless for a very 
limited range of the gradient of the temperature or when the gradients of the temperature 
of the electrons and ions have opposite signs which both cases are very unlikely. However, 
in the case when the heat flux does not play the role, the system can be unstable according 
to the Schwarzschild criterion. 

The current linear analysis is for simplified initial conditions, in which the background 
magnetic field, temperature gradient, and gravity are along the same direction. However, 
another configuration should also be examined using the £^-approach, in which the initial 
magnetic field is perpendicular to the direction of gravity. This will be done in the forth- 
coming paper. 
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